Coupling of bouncing-ball modes to the chaotic sea and their counting function 
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We study the coupling of bouncing-ball modes to chaotic modes in two-dimensional billiards with 
two parallel boundary segments. Analytically, we predict the corresponding decay rates using the 
fictitious integrable system approach. Agreement with numerically determined rates is found for 
the stadium and the cosine billiard. We use this result to predict the asymptotic behavior of the 
counting function Nbb(E) ~ E 5 . For the stadium billiard we find agreement with the previous result 
5 = 3/4. For the cosine billiard we derive 5 — 5/8, which is confirmed numerically and is well below 
the previously predicted upper bound S = 9/10. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

Two-dimensional billiard systems have found abun- 
dant applications in contemporary physics such as 
for electromagnetic and acoustic resonators, microdisk 
lasers, atomic matter waves in optical billiards, and quan- 
tum dots [iMSl. The classical dynamics of these billiards 
is described by a point particle of mass M moving with 
constant velocity inside a domain Q with elastic reflec- 
tions at its boundary dQ. Depending on the shape of the 
boundary the phase space can be regular, mixed regular- 
chaotic, or chaotic. Quantum mechanically, billiards are 
described by the time-independent Schrodinger equation 
(in units ft = 2M = 1) 



Aipi(x,y) = Eiijji{x,y), (x,y) G Q 



(1) 



with Dirichlet boundary condition ^ = on <90, eigen- 
functions -0/, eigenvalues Ei, Laplace operator A, and 

leN. 

Many nonintegrable billiards of interest contain a rect- 
angular region combined with other boundary segments, 
for example the stadium [6[ , the Sinai |7[ , and the cosine 
billiard |8[; see Fig. [TJ Classically, in these billiards a 
family of marginally stable periodic orbits exists, which 
bounce with perpendicular reflections between the two 
parallel parts of the boundary. These so-called bouncing- 
ball orbits are surrounded by chaotic motion in phase 
space. For simplicity we will only consider billiards which 
show no visible regular regions. Quantum mechanically, 
most of the eigenfunctions ipi are chaotic; i.e., they ex- 
tend over the whole chaotic phase space. In contrast, 
the so-called bouncing-ball modes ^ mn [9l-ll4j concen- 
trate on the marginally stable periodic orbits. They have 
a structure similar to the eigenstates of a rectangle; see 
Figs. QJc) andQJd), an d are characterized by the quan- 
tum numbers m and n, where m describes the quanti- 
zation along the bouncing-ball orbits and n perpendic- 
ular to them. In Ref. [15] it was proven that semiclas- 
sically the modes with n = 1 and increasing m exist. 
For the counting function N\y\y(E) of the bouncing-ball 
modes up to energy E this sequence of modes leads to 



A^bb(^) ~ \/~E. Typically also modes with higher exci- 
tations n > 1 perpendicular to the periodic orbits exist 
for large enough m. Which bouncing-ball modes ^ mn 
are realized for a specific billiard depends on the cou- 
plings of the bouncing-ball modes to the chaotic modes. 
For fixed m these couplings increase with n such that for 
large n the bouncing-ball modes couple to many neigh- 
boring chaotic modes and thus disappear. One expects 
that the number of bouncing-ball modes is asymptoti- 
cally described by a power law N\^{E) oc E 5 with ex- 
ponent 1/2 < S < 1. The exponent 6 = 1 cannot 
be achieved in ergodic billiards, as quantum ergodicity 
[l6l-[20j requires that the fraction of exceptional eigen- 
functions must vanish, i.e., Nhb(E)/N(E) — >• 0, where 
N(E) ~ E is the total number of eigenstates. 

The exponent S depends on the shape of the billiard 
in the vicinity of the rectangular bouncing-ball region. 
In Ref. [2l| it was shown that for the stadium billiard 
the exponent S = 3/4 arises, using an EBK-like quanti- 
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(c) 





FIG. 1. (Color online) Schematic pictures of the desym- 
metrized (a) stadium and (b) cosine billiard. Each billiard 
has a rectangular bouncing-ball region (gray shaded) of width 
w and height h. In the rectangular region so-called bouncing- 
ball orbits exist perpendicular to the parallel parts of the 
boundary (red vertical lines). In (c) and (d) bouncing-ball 
modes are shown for the stadium and the cosine billiard, re- 
spectively. 
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zation of the bouncing-ball modes. With a different ap- 
proach based on an adiabatic approximation [lOj of the 
bouncing-ball modes an upper bound for the exponent S 
was obtained for any chaotic billiard with a rectangular 
bouncing-ball region [22]. For the stadium billiard this 
bound agrees with (5 = 3/4 from Ref. [2l|. For the cosine 
billiard a bound of S = 9/10 was obtained and S « 0.87 
was observed numerically. 

In this paper we relate the couplings between 
bouncing-ball modes and chaotic modes to the num- 
ber Nbb(E) of bouncing-ball modes in a billiard. The 
couplings give rise to decay rates 7, which describe the 
initial exponential decay ~e _7t of states concentrating 
on the marginally stable periodic orbits to the chaotic 
region of phase space. In order to predict the decay 
rates 7 we employ the fictitious integrable system ap- 
proach [H, [24]], which was previously used to deter- 
mine regular-to-chaotic tunneling rates in systems with a 
mixed phase space [23l-l27| . We find a power-law decrease 
of the decay rates with increasing energy. We argue 
that the decreasing 7 imply the semiclassical existence 
of the bouncing-ball modes, complementing previous ap- 
proaches [111, LL5|, |2l|. Furthermore, we use this predic- 
tion of decay rates to count the number of bouncing-ball 
modes iVbb (-£")• For the stadium billiard we confirm the 
prediction of S = 3/4 [U, [22| . For the cosine billiard we 
derive the exponent 6 = 5/8 = 0.625, which is well be- 
low the previously predicted upper bound S = 9/10 [22j |. 
Numerically we observe 5 ~ 0.64 using a larger energy 
range than in Ref. [22|, where S « 0.87 was found. 

This paper is organized as follows. In Sec. [II] we review 
the fictitious integrable system approach, derive a pre- 
diction for the decay rates of bouncing-ball modes, and 
compare the results to numerically determined rates. As 
examples we consider the ergodic stadium billiard and 
the cosine billiard at parameters, where no regular mo- 
tion is visible in phase space. In Sec. [Till we use the pre- 
diction of decay rates to determine the number TVbb(^) 
of bouncing-ball modes and compare the results to the 
literature 0, [22|. A summary is given in Sec. IIVI 



II. COUPLING OF BOUNCING-BALL MODES 

For the analysis of the coupling of bouncing-ball modes 
to the chaotic sea we use the analogy to systems with a 
mixed phase space. There the coupling of regular and 
chaotic modes is caused by dynamical tunneling [28| , 
which quantum mechanically connects the classically sep- 
arated regular and chaotic regions. The billiards studied 
in this paper have no regular islands in phase space. How- 
ever, there is a region surrounding the marginally sta- 
ble bouncing-ball orbits which acts like a regular island, 
whose border is energy dependent such that its area semi- 
classically goes to zero [21[. For this situation we apply 
the fictitious integrable system approach [23], [24| to pre- 
dict the couplings of bouncing-ball modes to the chaotic 
sea. Whether this coupling occurs due to dynamical tun- 



neling as in mixed systems, or due to classically allowed 
transitions, e.g., through partial barriers, or some other 
mechanism is an open question that we do not address 
in this paper. 



A. The fictitious integrable system approach 

We first give a brief review of the fictitious integrable 
system approach [24| which was previously used to deter- 
mine regular-to-chaotic tunneling rates in systems with 
a mixed phase space [23l-[27|. Here it is applied in order 
to predict couplings of bouncing-ball modes in chaotic 
billiards. The main idea of the fictitious integrable 
system approach is the decomposition of Hilbert space 
into two parts, which correspond to the bouncing-ball 
modes (regular subspace) and to the remaining chaotic 
modes (chaotic subspace). Such a decomposition can be 
found by introducing a fictitious integrable system H reg . 
This system has to be chosen such that its classical dy- 
namics is integrable and contains the marginally stable 
bouncing-ball motion of the chaotic billiard. Quantum 
mechanically, the eigenstates IV^g ) °f ^reg? ^regl^g ) = 
E™™ l^reg)? closely resemble the bouncing-ball modes of 
the chaotic billiard and are used as a basis for the regular 
subspace. The regular modes IV^g ) are characterized by 
two quantum numbers m and n, where m describes the 
quantization along the bouncing-ball orbits and n per- 
pendicular to them. They localize on quantizing tori of 
i^ reg , and decay beyond. This decay is the decisive prop- 
erty of |^reg )> which have no chaotic admixture, in con- 
trast to bouncing-ball modes in a chaotic billiard. Two 
choices for the system H reg will be discussed in Sec. Ill Bt 

Introducing a basis l^ch) m the chaotic subspace, the 
coupling matrix element of one regular mode and 
a chaotic mode |^ch) is given by 



(2) 



The disadvantage of these couplings v c h,mn is their de- 
pendence on the size of the chaotic region via the nor- 
malization of the chaotic modes |^ch)- A better suited 
quantity, that is not affected by such changes, is the decay 
rate 7. Under variation of the billiard boundary far away 
from the rectangular bouncing-ball region it is unique for 
each regular state From the couplings v c h, mn the 

decay rate is obtained with Fermi's golden rule 



7mn = 27r(|v c h, mn | 2 )Pch- 



(3) 



Here we average over the modulus squared of the fluctu- 
ating coupling matrix elements v c h,mn of one particular 
regular mode and different chaotic modes of similar en- 
ergy. The chaotic density of states is approximated by 
the leading Weyl term p c h ~ A/ (47r), in which A denotes 
the area of the billiard. 

To illustrate the decay rates 7 one may consider a reg- 
ular state l^reg ) concentrating close to the marginally 
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FIG. 2. (Color online) Comparison of bouncing-ball modes of (a) the cosine billiard with quantum number (23, 4) and (b) the 
stadium billiard with quantum number (23, 5) to the corresponding adiabatic modes [(c), (d)] and to the rectangle modes [(e), 
(f)]. The adiabatic modes (c) and (d) closely resemble the bouncing-ball modes (a) and (b), respectively. The rectangle mode 
(e) closely resembles the bouncing-ball mode (a) of the cosine billiard, while for the stadium billiard the rectangle mode (f) 
fails to reproduce the bouncing-ball mode (b). 



stable bouncing-ball orbits which is coupled to a contin- 
uum of chaotic states. Its decay ~ e -owt is character- 
ized by the rate 7 mn . For systems with a finite phase 
space this exponential decay occurs at most up to the 
Heisenberg time th = 27r/A, where A is the mean level 
spacing. Alternatively, the decay rates are the inverse of 
the lifetimes of resonances in a corresponding open sys- 
tem, which can be obtained, e.g., by adding an absorbing 
region in the chaotic component of phase space far away 
from the bouncing-ball region. 



B. Integrable approximations 

In the following we predict decay rates 7 of bouncing- 
ball modes with Eqs. (j2j) and (j3j). For this we con- 
struct a fictitious integrable system i^ reg , whose eigen- 
states |^^g ) resemble the bouncing-ball modes of the 
chaotic billiard. Two approaches will be presented be- 
low. The first uses adiabatic modes, which approximate 
the bouncing-ball modes quite well, but have to be eval- 
uated numerically. The second approach uses modes of 
the rectangular billiard, which are given analytically, but 
turn out to be a good approximation for the bouncing- 
ball modes of the cosine billiard only. 



1. Adiabatic-mode approximation 

In order to construct modes which approximate the 
bouncing-ball modes of a billiard we use the adiabatic 



separation ansatz [Io| 
with 



2 . / m7ry\ 

sm [Wj ) 



h(x) 



(4) 



(5) 



where h(x) is the height of the billiard at x with h(0) = h. 
Equation (j4j is close to a separation ansatz but (f m 
weakly depends on x in order to satisfy the Dirichlet 
boundary condition on dfl. The function cp m accounts for 
the quantization of the fast bouncing-ball motion in the 
y direction. The slow motion in the x direction is quan- 
tized by demanding that Xmn fulfills the ID Schrodinger 
equation [10[, 



iXmn(x), (6) 



which we solve numerically. Here the effective potential 
is m 2 V(x) with 



V(x) = 



h{xf 



7T 



(7) 



and the eigenvalues are denoted by e mn . From c mn 
one obtains approximate eigenvalues of the bouncing-ball 
modes E™™ = m 2 7r 2 /h 2 + e mn . The ansatz, Eq. (|4]), ex- 
actly satisfies the Schrodinger equation inside the rect- 
angular region of the billiard. Outside of the rectangular 
region this is only approximately true. There the func- 
tions Xmn(x) will eventually decay exponentially, as the 
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effective potential m 2 V(x) becomes very steep for large 
m. 

The adiabatic modes show good agreement to the 
bouncing-ball modes of the stadium and the cosine bil- 
liard if n <C m; see Figs. Efa)-|2fd). In particular the 
bouncing-ball modes are well reproduced in the region 
x ~ w, where the rectangular region ends. Note that 
small deviations arise for modes with large n, in particu- 
lar for the stadium billiard (not shown). However, these 
modes strongly couple to chaotic modes and thus do not 
contribute to the counting function in Sec. IIIII 

We now predict decay rates of bouncing-ball modes 
using the adiabatic modes l^^g )• Replacing the coupling 
matrix elements v c h,mn in Eq. (j3j) by Eq. (j2j) we find for 
the decay rates 



27rp ch 
N A 

27rp c h 



<vc g n i# ( E i^ch)(^hi]^iCg n )(8) 

ch,AE 



\PcKAeH\^)\\ 2 - 



(9) 



To obtain Eq. flHJ) we express the average in Eq. (|2j) by a 
sum over all chaotic states |^ch) in an energy interval AE 
around the energy E™™ of the adiabatic mode and Na is 
the number of chaotic states in this energy interval. The 
term ^2 ch AE |V ; ch)(V ; ch| in Eq. flSJ is equal to a projector 
P c h,AE onto the chaotic subspace in the considered en- 
ergy interval. We approximate this projector using the 
adiabatic modes |^™g), 



PchAE « Pae 1 



E 

m'>n' 



(10) 



PaeH\i/>™) = Ea E ciEi\^), where c, = 
With this procedure Eq. ([TT]) can be evaluated numer- 
ically for all bouncing-ball modes of quantum number 
(m,n). We find good agreement to numerically deter- 
mined rates for the cosine and the stadium billiard; see 
Sec. ECl 



2. Rectangle-mode approximation 

While the approach for the determination of decay 
rates presented in the last section is applicable for any 
billiard with a rectangular bouncing-ball region, an an- 
alytical result is not available. In order to find such an 
analytical result we now approximate the bouncing-ball 
modes by the eigenstates l^g ) °f a rectangular billiard 
and construct the chaotic states |^ch) entering in Eq. (J2j) 
by a random- wave model [29]. As chaotic states l^ch) 
modeled in this way are not orthogonal to the regular 
states l^^), Eq. ([2]) has to be modified (see Sec. II A in 
Ref. [24[), which leads to the approximation 



^Ch ' 



(^ch|# - ^regl^i 



reg / 



(13) 



The simplest approximation of a rectangular bouncing- 
ball region of width w and height h is given by the rect- 
angular billiard of the same width and height. We now 
use this billiard as the fictitious integrable system H reg 
in Eq. ([T3]) . Its eigenstates in position representation are 
given analytically by 



:Sin l— J sin l 



miry\ 



(14) 



-pmn 
J reg 



Here Pae projects onto the energy interval [E" ( 
AE/2,E™ + AE/2}. The choice m' > n' for the sum 
ensures that P c h,A£; projects onto the phase-space com- 
ponent with k x /h > ky/w, which excludes the bouncing- 
ball region but also parts of the surrounding chaotic re- 
gion. This choice will underestimate the decay rates 
slightly. However, a more precise projection would re- 
quire an a priori knowledge of the energy-dependent bor- 
der between the bouncing-ball region and the chaotic re- 
gion. Using Eq. ([TP]) in Eq. (J9j) the decay rates are de- 
termined by 



with eigenenergies E\ 



2Na 



PaeH\^}~ E ^W,mn|Cg U/ ) 
(m' ,n')eAE 



where 



v. 



m'n' ,mn 



reg / 



(ii) 



(12) 



and p c h ~ A/(4tt) has been used. In order to numeri- 
cally evaluate the term PaeH^ 7 ^) in Eq. ([TT]) we ex- 
pand the state |^reg ) m the eigenbasis |^) of H in the 
energy interval [E™ - AE/2,E™ + AE/2] such that 



ran 
reg 



7T 2 m 2 /h 2 + 7T 2 n 2 /w 2 . Note 



that these rectangle modes ip™™(x, y) are zero for x > 
while the bouncing-ball modes of chaotic billiards will 
typically extend into this region. Hence, the rectangle 
modes will only be an approximation of the true eigen- 
states of the billiard. Figure [2] shows that the rectan- 
gle modes are good approximations for the bouncing-ball 
modes of the cosine billiard, which has a corner of angle 
7r/2 at x = w, but not for the bouncing-ball modes of 
the stadium billiard, which substantially extend beyond 
x = w. As the fictitious integrable system approach re- 
lies on a good approximation of the bouncing-ball modes, 
the following prediction of decay rates is expected to ap- 
ply only for billiards which have a corner of angle ir/2 at 
x = w, like the cosine billiard. 

When evaluating Eq. ([T3]) an infinite potential differ- 
ence arises between the Hamiltonian H of the chaotic 
billiard and H reg of the rectangle, H — H reg — — oo, 
for x > w. At the same time for the rectangle modes 
i/j™g(x,y) = holds in that region, which leads to an 
undefined product — oo • 0. Similar to the approach pre- 
sented for the mushroom billiard in Refs. [24], |25| we cir- 
cumvent this problem by introducing a rectangular bil- 
liard which is extended by a finite potential W at x > w 
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for which in the end the limit W —> oo is taken. This 
leads to 



dyipch(x w,y)d x ^™(x = w,y). (15) 



For completeness we now give the derivation of Eq. ([15]) 
following Refs. [24], [25j: We first introduce a rectangular 
billiard which is extended by a finite potential W for 

x > w, 

HW g =p 2 x +p 2 y + V(x : y) (16) 

{0 for < x < w, < y < h 
W for x > w, < y < h (17) 
oo otherwise 

and consider the limit W — » oo in which the original 
rectangular billiard is recovered. For finite W the regular 
eigenfunctions ^™g^w( x iy) °f -^reg decay into the region 
x > w, which is described by 

^Iwi^v) = ^Zwi* = w,y)e- x ^- w K (18) 

Here A depends on W via the Schrodinger equation as 

(19) 



A 2 + W = E w 



Since the derivative of the regular eigenfunctions 
has to be continuous at x = w we obtain 



ran 
reg,Vl^ 



Wg>(* = w,y) = -A ^™>(i = w, y). (20) 



This can be used for rewriting Eq. (JT3J) for the coupling 
matrix elements 



w +b h(x) 



= lim 

W^oo 

w 

The term 



lirn^y dx J dy^(x,y)(-W)^ iW {x,y)(21) 

w 
w+b h(x) 

Jdxjdy i, ch (x, y)^e- x ^d x ^ w (w, y).(22) 



W X (x-w) 

A 



W 



-^/W-E^ n {x-w) ^3) 



reduces in the limit W —> oo to 25 (x — w) for x > 
where we used Eq. (p~9|) and that E^ n remains bounded. 
This finally leads to Eq. (PS). 

Evaluating the coupling matrix elements, Eq. (fT5]h for 
the rectangle modes, Eq. JEJ), gives 



= (-ir- 



2 titt 



o\yip ch (x = sin — — . 



V h J 



Therefore the decay rates 7 mn , Eq. (|3)), scale with 



(24) 



In order to find the scaling of the decay rates with 
m we employ a random-wave model [29| for the chaotic 
states ip c h(x,y) entering in Eq. ([24]) . Here it is given by 

2 s max 571" \ 

ipch(x, y) = —j= ^ c s sin (—J sin(k XtS x + <p s ) (25) 



and accounts for the Dirichlet boundary conditions at 
y = and y = h, but ignores the corner at x = w. Here 
the c s are Gaussian-distributed random variables with 
mean zero and variance <j s = {c 2 s ) . The phases cp s are uni- 
formly distributed in [0,2tt), fc^ = y^g" ~7T 2 s 2 /h 2 , 

and 5 max = L(^/ 7r )v // ^reg'J • Furthermore we define 
k VyS = irs/h and the angles a s = arctan(/c 2/)S /A: ;E)S ). 
These angles are not equidistributed on [0,7r/2]; see 
Fig. El To approximate an equidistribution of direc- 
tions we compensate this by appropriate weights a s , with 
Ss=i x cr s = 1. A natural choice is given by assigning to 
each <j s half of the size of the adjacent angular intervals 



1 <C S <C 5 max 

-i)/2 , s — s max 

8 = 1, 



l f (ats+i - a s -i)/2 , 

n / 2 { (ai + a 2 )/2, 

(26) 

see Fig. [3] 

Using Eq. ([25]) in Eq. ([M]) the ?/ integration over the 
two sine functions gives S sm h/2 such that 

2 TITT h 

Vc^mn = {-l) n —^= -=c m sm(k x , m w + <p m ). (27) 



We now determine the decay rates with Eq. ([3]) using an 
ensemble average of the modulus squared of the coupling 




FIG. 3. (Color online) Illustration of the discrete directions 
contributing to the random- wave model, Eq. ([25]) . in k space 
at energy E™™ (black quarter circle) with m — 10 and n — 
2. While the wave numbers k y , s are equispaced, the wave 
numbers k x , s as well as the angles a s are not equispaced. 
The angular regions around each a s are marked by gray lines 
and the one for a SrTiaY is shaded. 
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matrix elements given by Eq. ([27]) . This results in 



In 



7T 2 h 



W° 



n 2 a n 



(28) 



which contains the weight <r m . In order to find the weight 
cr m from Eq. ([26]) we have to evaluate s max . One finds 
Smax = m f° r small enough n < \/2(m + 1/2) Z /ft fol- 
lowing from (h/7r)y/E™™ < m + 1. We now perform 
a first-order expansion of the angles a m in Eq. ([26]) for 
large m and restrict to even smaller n <C \/m. We obtain 



V2 1 



leading to our final result 



7r7 



V^Trft n 2 



(29) 



(30) 



which predicts the decay rates of bouncing-ball modes 
characterized by the quantum numbers m and n to the 
chaotic sea for n <C y/m. For increasing m the decay 
rates decrease like a square root while for increasing n 
they increase quadratic ally. 

The power-law decrease of the decay rates for constant 
n and increasing m is in contrast to the exponential de- 
crease found for direct regular-to-chaotic tunneling rates 
in mixed systems, as in the mushroom billiard [25], [30| 
and in quantum maps [23|, I3ll - l33| . The decreasing decay 
rates for fixed n lead to decreasing couplings v c h,mn, see 
Eq. (j3j), between the bouncing-ball modes and chaotic 
modes compared to the mean level spacing A. This im- 
plies for large enough m the semiclassical existence of 
the bouncing-ball mode (m, n). Note that the prediction, 
Eq. ([30]) . only depends on the rectangular bouncing-ball 
region but not on the other regions of the considered bil- 
liard. Hence, it is the same for each billiard studied in 
this paper. However, it can only be applied if the rect- 
angle modes closely resemble the bouncing-ball modes. 
This is the case for the cosine billiard but not for the sta- 
dium billiard; see Sec. Ill CI for the comparison of Eq. ([30]) 
to numerically determined rates. 



C. Applications 

We now evaluate Eqs. ([TT]) and ([30]) for the decay 
rates of bouncing-ball modes, as discussed in Sec. IIIB[ 
and compare the results to numerically determined decay 
rates for specific example systems. We choose the desym- 
metrized cosine billiard with w = 1, ft = 0.6, b = 0.2, 
h(x) = ftarccos(2(x — w)/b — 1)/tt for x > w, and la- 
bel this billiard by B c . Furthermore we consider the 
desymme trized stadium billiard with w = 1, ft = 0.6, 
h(x) = y/h 2 — (x — w) 2 for x > w, and label this billiard 
by B a . 

For the evaluation of Eq. ([TT]) the energy interval AE 
has to be specified. We choose AE w 100A with mean 
level spacing A. Note that for small energy intervals 



AE < 20A the results of Eq. ([TT]) show strong fluctu- 
ations, as the adiabatic states cannot be properly ex- 
panded in the eigenbasis of H. For AE > 20 A, however, 
the results show only small changes. 

In order to calculate decay rates numerically we de- 
termine the spectrum of the stadium and the cosine bil- 
liard, using the improved method of particular solutions 
[34j , under variation of parts of the billiard boundary: 
For the cosine billiard we vary the width of the cosine 
part b around b = 0.2 and for the stadium we vary w 
around w = 1. For small variations of these parame- 
ters the eigenvalues of the bouncing-ball modes remain 
almost unaffected while the eigenvalues of the chaotic 
states show strong variations, due to the changing den- 
sity of chaotic states. Analyzing up to 30 avoided cross- 
ings A£? c h, mn = 2v c h, mn of a given bouncing-ball mode 
Tprnn with chaotic states we deduce the decay rate from 
Fermi's golden rule pil |25|, 



7r7 



= l(\AE ( 



ch.mn I 



2 A). 



(31) 



[7] 



Here we average over all numerically determined widths 
AE c \ l ^ rnn and the corresponding billiard areas A. 

In Fig. |4ja) we show for the cosine billiard B c numeri- 
cal rates (dots) compared with the prediction of Eq. ([TT]) . 
using the adiabatic modes (solid lines), and Eq. ([30]) . 
using the rectangle modes (dashed lines). We consider 
bouncing-ball modes with horizontal quantum number 
n = 1,2 and increasing vertical quantum number m. 
We find very good agreement with deviations smaller 
than a factor of two for both predictions. The numer- 
ical rates and the results of Eq. ([TT]) follow the power 
law 7 ~ m -1 / 2 predicted by Eq. ([30]) . This confirms 
that for the cosine billiard the adiabatic modes as well 
as the rectangle modes sufficiently well approximate the 
bouncing-ball modes, as seen in Figs.[2ja),[2fc), and[2je 
Note that also for the desymmetrized Sinai billiard 
both predictions agree with numerical rates (not shown) . 
It also has a corner of angle tt/2 at x = w, such that 
the rectangle and the adiabatic modes closely resemble 
its bouncing-ball modes. 

Figure SJb) shows the decay rates for the stadium bil- 
liard B s for n = 1, 2 and increasing m. We compare the 
numerical rates (dots) to the prediction of Eq. ([TT]) (solid 
lines), and Eq. ([30]) (dashed lines). Here, the average be- 
havior of the fluctuating numerical rates agrees with the 
prediction of Eq. ([TT]) . As for the cosine billiard this pre- 
diction seems to follow a power-law for large m, however, 
with an exponent smaller than —1/2, close to —1. As 
expected, Eq. ([30]) . based on the rectangle modes, does 
not reproduce the numerical rates, showing deviations 
by a factor of 10. These deviations arise as the rectangle 
modes do not agree well enough with the bouncing-ball 
modes of the stadium billiard, in contrast to the adia- 
batic modes used in Eq. ([TT]) ; see Figs. [2fb), [2fd), and 
[2ff). Note that the small fluctuations in the prediction 
arise due to the numerical evaluation of Eq. ([TT]) and 
might be related to the chosen projector P c h,AE- The 
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FIG. 4. (Color online) Decay rates of bouncing-ball modes 
for (a) the cosine billiard B c and (b) the stadium billiard 
B s . For the quantum numbers n = 1,2 and increasing m 
we compare numerical rates (dots) to the results of Eqs. (|11[) 
(solid lines) and (|30|) (dashed lines) on a double-logarithmic 
scale. The insets in (a) and (b) show the bouncing-ball modes 
of quantum number (7, 1) and (39, 2). 



large fluctuations in the numerically determined rates 
could be caused by additional couplings of the bouncing- 
ball modes to scars concentrating on unstable periodic 
orbits, which is left for future studies. 

In Fig. Owe show numerically determined decay rates 
and the two predictions, Eqs. ([TT]) and ([30]) . for bouncing- 
ball modes with fixed quantum number m = 30 and in- 
creasing n. For both the cosine and the stadium billiard 
we find good agreement between the numerical rates and 
Eq. ([TT]) using the adiabatic modes. Equation ([30]) , based 
on the rectangle modes, is valid only for the cosine bil- 
liard as discussed before. For increasing n the decay rates 
increase almost quadratically, 7 ~ n 2 , for both of the bil- 
liards, as found in Eq. ([30]) . 



(a) (b) 




FIG. 5. (Color online) Decay rates of bouncing-ball modes for 
(a) the cosine billiard B c and (b) the stadium billiard B s . For 
the quantum numbers rn = 30 and n — 1,2,..., we compare 
numerical rates (dots) to the results of Eqs. (|lip (solid lines) 
and ([30]) (dashed lines) on a double-logarithmic scale. The 
insets show the bouncing-ball modes of quantum number (a) 
(30,4) and (b) (30,6). 



III. COUNTING BOUNCING-BALL MODES 

We now study and predict the number of bouncing- 
ball modes N\^{E) up to energy E for a billiard with a 
rectangular region. Due to the quantum ergodicity the- 
orem [16] in an ergodic billiard the ratio of N^(E) to 
the total number of eigenstates N(E) goes to zero in the 
semiclassical limit, N\ ) \ ) (E)/N(E) — >• 0. In addition the 
number of bouncing-ball modes should increase at least 
as V^E, which gives the number of states concentrating 
on a one-dimensional line. Hence, one expects for the 
number of bouncing-ball modes 

N hh (E) = a-E 5 + ... (32) 

with exponent 1/2 < 5 < 1 and prefactor a. Here we 
neglect higher order terms which contain powers of the 
energy smaller than S. 

Note that in Ref. [35j the boundary contribution ~\^E 
of Eq. ([32]) has been determined for subsets of eigenstates 
in billiards. For the bouncing-ball modes in the billiards 
considered in this paper this boundary contribution sug- 
gests the additional term —h/(Air)\fE. As we are not 
dealing with a regular region around a fixed point but an 
energy dependent region around the line of marginally 
stable bouncing-ball orbits, we expect the boundary con- 
tribution — 2/i/ (47r) \[E. However, it will not be consid- 
ered in the following, as we are interested in the leading 
term only. 



A. Counting using overlap with adiabatic modes 

First we count the number of bouncing-ball modes 
N\y\)(E) in the cosine and the stadium billiard numer- 
ically. For this purpose we calculate their first ~3000 
eigenstates and eigenvalues E\ using the improved 
method of particular solutions [34J. We then determine 
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the overlap 

w = mm^)\ 2 (33) 

of each eigenstate with the adiabatic approximations 
l^reg ) °f the bouncing-ball modes. If |^) has an overlap 
W > 1/2 with an adiabatic mode with m > n we consider 
l^z) as a bouncing-ball mode of the system. Counting 
the number of these modes gives a numerical estimate of 
N hh (E). 

In Figs.[6fa) and[6jb) the results (black staircase func- 
tions) for the cosine and the stadium billiard are shown, 
respectively, on a double-logarithmic scale. For large 
energies the data follow in both cases a power law for 
the number of bouncing-ball modes, as expected from 
Eq. ([32]). Note that varying the cut off W = 1/2 in 
[0.5,0.6] changes the prefactors a but not the exponents 
5 (not shown). 

Specific ally for the stadium billiard one expects S = 
3/4 HH, HI- In Fig. [6{b) we show the function 
N hh (E) = 0.13 • £ 3/4 (red dashed line), which on a 
double-logarithmic scale has the same slope as the nu- 
merical data. In the inset we fit an exponent 5 ~ 0.73 to 
the data which is close to S = 3/4. 

For the cosine billiard S — 9/10 was obtained as an 
upper bound and S ~ 0.87 was observed numerically 
[22[ | . However, in Fig. [6fa) the data are best fitted by 
an exponent 5 « 0.64 (red dashed line in the inset). In 
order to resolve this contradiction we also studied the 
cosine billiard used in Ref. [22| with parameters h = 1, 
w = 2, and b = 1 (not shown). In Ref. [22| numerical 
data for N^^E) is presented which is well described by 
A/"bb(£0 = &E 5 + /3 with an exponent S « 0.87 and a 
constant offset /?. For the identification of the bouncing- 
ball modes a visual selection is performed and states up 
to energy E = 8400 are considered. If we use the over- 
lap criterion with W > 1/2 for this cosine billiard and 
consider states up to energy E = 20000 we obtain an ex- 
ponent S ~ 0.64. Also a visual selection over this larger 
energy range gives an exponent S « 0.7. Hence, we con- 
clude that the smaller energy interval used in Ref. [22| is 
not sufficient to find the asymptotic scaling exponent. 



(a) 
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FIG. 6. (Color online) Number of bouncing-ball modes 
Nhb(E) for (a) the cosine billiard B c and (b) the stadium 
billiard B s > We compare the results using the overlap crite- 
rion (black staircase functions) and the prediction using the 
decay rates (gray staircase functions) to power laws with ex- 
ponent (a) 6 = 5/8 and (b) 6 = 3/4 (red dashed lines) on a 
double-logarithmic scale. In (a) we show in addition a power 
law with exponent S = 9/10 (dotted line) [22j]. The insets 
show the results using the overlap criterion on a linear scale 
compared to Eq. (j32|) fitted for E > 10000 (red dashed lines). 
We find S = 0.64, a = 0.196 in (a) and S = 0.73, a = 0.112 
in (b). 



B. Counting using decay rates 

In Refs. [36|, [37j the criterion 7 mn < 1/th has been 
found to describe the existence of regular states for sys- 
tems with a mixed phase space. It compares the decay 
rates 7 mn to the Heisenberg time th = h/A = 2tt/A. 
Here we apply a similar criterion to predict the number 
of bouncing-ball modes in a billiard, using the results 
([HD and ([30]) for the decay rates derived in Sec. [TTJ 

From Eq. (|3|) we obtain the average coupling v mn of a 
bouncing-ball mode to the chaotic sea, 

Vmn = V w (34) 



If this coupling v mn is much smaller than A, on average 
the bouncing-ball mode weakly couples to the chaotic 
modes such that it is visible. If the coupling is much 
larger than A, on average the bouncing-ball mode couples 
strongly to many chaotic modes such that it is not visi- 
ble. In between these two limiting cases there is a smooth 
transition, which was studied in Ref. [37[ for systems with 
a mixed phase space. In order to count the number of 
bouncing-ball modes we approximate this smooth tran- 
sition by a sharp condition. The criterion 7 mn < 1/th 
would lead to the condition v mn < A/(2tt). It is very 
strict in the sense that it only allows for very small cou- 
plings between the bouncing-ball modes and the chaotic 
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modes. Approximately in the middle of the transition we 
find the condition 

Vmn < ^ ( 35 ) 

which we use in the following for counting the number of 
bouncing-ball modes. 

We now calculate the average coupling v mn for all 
bouncing-ball modes up to energy E with Eq. ([34]) , using 
the decay rates j mn from Eq. ([TT]) . We then count the 
number of those modes which fulfill the condition ([35]) . 
The results are shown in Fig. [6ja) for the cosine billiard 
and in Fig. [6^b) for the stadium billiard as gray stair- 
case functions. For the stadium billiard we obtain the 
exponent S « 0.76 close to S = 3/4 [2l|, [22| and for the 
cosine billiard we find S « 0.63. Note that changing the 
condition ([35]) to v mn < A/(2ty) changes the prefactor a 
but not the exponent 5 for the examples we considered 
(not shown). 

In order to obtain an analytical prediction of N^(E) 
for the cosine billiard we use the result for the decay rates 
7 mn , Eq. ([30]). Together with Eq. ([34]), condition ([35]) is 
fulfilled by all bouncing-ball modes with 

/ 7rA W 3 i , . 

n< y-iv^h ml - (36) 

For large m this complies with the restriction n <C \pm 
used in the derivation of Eq. ([30]) . For the determination 
of N\fo(E) one now has to count all states of quantum 
number (m, n) for which Eq. ([36]) is fulfilled and E™™ = 
7r 2 m 2 /h 2 +7r 2 n 2 /w 2 < E. This sum can be approximated 
by an integral, where we integrate the right-hand side of 
Eq. ([36]) over m in the interval [0, ra max ] and approximate 
m max « yfEhj-K. This finally gives for the number of 
bouncing-ball modes 

N hh (E)^a-E~s. (37) 

with the constant a = 4w 3 ^ 2 h 3 ^y^A/2/(5 • 2 1 / 4 tt 7 / 4 ). 
Equation ([37]) gives the exponent J = 5/8 [red dashed line 
in Fig. |B(a)]. It is close to the exponents (5 « 0.64 and 



(5 ~ 0.63 which were obtained using the overlap criterion 
and the decay rates from Eq. ([TT]) , respectively. Even the 
prefactor a « 0.34 agrees roughly with the fitted value 
a ~ 0.2. Note that Eq. ([37]) is only valid for billiards 
with a confining corner at x = w, as in the cosine billiard, 
because otherwise Eq. ([30]) for the decay rates cannot be 
applied, as in the case of the stadium billiard. 

IV. SUMMARY 

In this paper we study the couplings of bouncing-ball 
modes to chaotic modes in two-dimensional billiards and 
use these couplings to count the number of bouncing-ball 
modes. In Sec. ITT1 we apply the fictitious integrable sys- 
tem approach [23|, |24[ to predict decay rates 7 mn , which 
describe the initial decay of bouncing-ball modes into the 
chaotic sea. Using the adiabatic modes as approximate 
bouncing-ball modes we evaluate Eq. ([TT]) and find agree- 
ment to numerical rates for the cosine and the stadium 
billiard. For the cosine billiard, which has a corner of 
angle tt/2 at x = we evaluate the fictitious integrable 
system approach analytically using the eigenstates of the 
rectangular billiard as approximate bouncing-ball modes. 
This leads to Eq. ([30]) which shows excellent agreement 
with numerical rates. As a result we find that the decay 
rates of bouncing-ball modes of constant quantum num- 
ber n decrease as a power law with increasing energy. 

In Sec. IIIII we use the results on the decay rates 
in order to count the number of bouncing-ball modes 
N hh (E) ~ E 5 up to energy E. For the stadium billiard 
we recover the exponent (5 = 3/4 [2l|, [22|. For the co- 
sine billiard we find 5 ~ 0.64. Using the analytical result 
Eq. ([30]) for the decay rates we derive 6 = 5/8, which 
is in agreement with our numerics and well below the 
previously predicted upper bound 6 = 9/10. 
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